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We consider a model of a thin film elastically attached to a rigid substrate. In the case in which 
the film expands relative to the substrate and assuming certain non-linear elastic behavior of the 
film, expansion ridges may appear, in which the material has collapsed, and the density is higher. By 
studying numerically this process, the possible morphologies of these collapsed regions are presented. 
They range from circular spots and straight stripes, to wiggle polygonal patterns and ring-shaped 
domains. The similarity of some of these results with patterns observed in delamination of thin 
films and bi-phase epitaxial growth is emphasized. 



PACS numbers: 

I. INTRODUCTION 

When the surface of a material contracts with respect 
to the underlying part, tensile stresses appear in it that 
can lead to the formation of a crack pattern onto the 
surface 0. Mud cracking can be considered to be the pro- 
totype of this surface fragmentation process. The main 
ingredients of surface fragmentation are a (quasi) two- 
dimensional film attached to a substrate (we always as- 
sume here the substrate is rigid) and a greater expansion 
of the substrate compared to the film as a function of 
some external control variable (typically humidity con- 
centration or temperature). The film and substrate can 
be the same material, as in mud cracks, and in that case 
it is only a difference in humidity concentration or tem- 
perature what identifies film and substrate. 

We will concentrate here in a process that in a certain 
sense is the inverse of surface fragmentation: we consider 
the film expanding with respect to the substrate. Due 
to the coupling to the substrate, the greater expansion 
of the film generates compressive stresses into it. No im- 
portant effects are expected if the film responds linear 
elastically to any deformation. But if the film can col- 
lapse upon compression, a coexistence of collapsed and 
non-collapsed regions in the film is expected. The pre- 
sentation and discussion of the different morphologies of 
these collapsed regions are the main aim of this paper. 

The two main features that fix the morphology of col- 
lapsed regions are the degree of mismatch between film 
and substrate, and the kind of non-linear elastic behav- 
ior of the film, in particular, the characteristics of the 
collapsed state. Different possibilities are studied here. 
Most of the results to be presented correspond to the 
case in which the system is isotropic in the plane of the 
film (the x-y plane), but some results for a model with 
square symmetry will also be presented. In the next Sec- 
tion we present the model and details on the simulation 
technique. In Section III we present the main results, 
and in Section IV we discuss two experimental situations 
in which the present model can be applied, namely, de- 
lamination patterns of thin film and structures appear- 
ing during bi-phase epitaxial growth. Section V contains 
some summary and conclusions. 



II. THE MODEL 

The model to be used is an extension of that used in 
0,0 to describe fracture (see also Q and 0)- It consid- 
ers the film in a two-dimensional approximation as de- 
scribed by the horizontal displacement field u(r). From 
these variables the two dimensional strain tensor is calcu- 
lated as Eij = l/2(dui/dxj + duj/dxi). For convenience, 
instead of we will use the following variables 

ei = (en + £22)/2 

e 2 = (en-e a2 )/2 (1) 

£3 = £12 = £21 

These three variables are not independent. They satisfy 
the St. Venant compatibility constraint 0, 0, |g, which 
reads 

{81 + dl)e x - {81 - d 2 y )e 2 - 2d x d y e 3 = 0. (2) 

The free energy density of the system contains three 
terms: a local term /o, a gradient term /y, and a sub- 
strate interaction term f su b s - The existence of a col- 
lapse transition for the film is encoded in the form of 
the local free energy term /q. If we intend to describe 
an isotropic material, only rotationally invariant combi- 
nations of the basic variables e, should enter the free 
energy. Those that can be constructed from the ej are 
ei itself, and = \/e§ + e|. A perfectly elastic mate- 
rial is described by a local free energy having a quadratic 
minimum, namely its local free energy is of the form 

ft Me = B(ei-4) 2 +n{4 + 4) (3) 

where B and [i are respectively proportional to the bulk 
and shear modulus of the material. With the present 
choice, in a completely relaxed state the system has e\ — 
e?, e 2 = e 3 = 0. 

If the material can collapse, fo must have some other 
minimum at some ef , describing the collapsed state. 
The position of the collapsed minimum (in particular, the 
ratio v = e^/(ef — e")) will have important consequences 
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in the morphologies of the collapsed regions that will be 
observed. Different locations of the collapsed minimum 
will be explored and discussed here. They are qualita- 
tively depicted in Fig. as cases A, B, C, and D. In 
most cases the morphologies of the collapsed patterns 
depend only on the position of the collapsed minimum 
with respect to the elastic minimum, the detailed form 
of fo being of minor importance. Only in case A, there 
are important variations depending on the value of the 
shear modulus of the collapsed state pF with respect to 
the value in the normal phase /x. Thus, two cases will 
be distinguished: Al, in which fi c is greater or equal p, 
and A2, in which pF is lower than p. Note that when 
the collapsed state has ^ 0, it corresponds actually to 
a ring of minima in the ei, e2, e-$ space. In other words, 
the collapse of a circular piece of material to a state with 
ed 7^ produces an elongated object, but the orientation 
of this object in the x-y plane can be any (in the draw- 
ings of Fig. this direction was chosen as the vertical 
one). The explicit form of the /o part of the free energy 
for the cases A-to-D is given in the appendix. 

Gradient terms in the free energy will be taken in the 
form 

/ v = Yl «*( Ve *) 2 > ( 4 ) 

i=l,2,3 

where ct2 = ot^ should be chosen to retain rotational in- 
variance. In all results below we take all Oj = a. 

The elastic interaction with the substrate is easily writ- 
ten in terms of the displacement variables u: 

f subs = l\u(r)\ 2 (5) 

where 7 measures the stiffness of this interaction. As we 
take the components of e to be our basic variables, we 
have to recast this energy in terms of them. This is more 
easily done in the Fourier space, and the result can be 
easily written after integration over the whole system as 

Where eij(k) are the Fourier transforms of the original 
ej(r). Here the tilde in the integral indicates that the 
k = mode is excluded. To avoid a divergent energy 
contribution, the value of ei(k = = 0) = e± (where the bar 
notes a spatial average) has to be adjusted for the system 
to fit on average to the substrate, namely, the mismatch 
with the substrate is incorporated in the model precisely 
through the value given to e%. 

The existence of non trivial spatial patterns in the 
model appears as a consequence of a competition between 
gradient and substrate term. In fact, in Fourier space a 
mode of wave vector ~ k and amplitude 1 produces a 
contribution to the energy of the order of ak 2 from the 
gradient terms, and "f/k 2 from the substrate term. The 
sum of this two contribution has a minimum at a value of 
k of the order of (7/a) 1 / 4 . This is the order of magnitude 



of the main spatial variations that will be seen to appear 
in the simulations. 

The equations of motion are taken to be of the over- 
damped form, namely 

de l SF , 

~~dT = <fei < <=1 ' 2 ' 3 > (?) 

where 

F = J dr 2 (f + / v + f subs ) . (8) 

The Saint Venant constraint is implemented by the use 
of Lagrange multipliers. Full details can be seen in 2] 
and H- 

III. RESULTS 

For each set of parameters, the system will settle down 
in a configuration that minimizes the total energy of the 
system. We will see that in general, metastable states ap- 
pear very often in the simulations. In order to get as close 
as possible to the true ground state configuration, an an- 
nealing process was implemented in which a stochastic 
term was added to the right hand side of Eq. Q, and 
the intensity of this term was progressively reduced down 
to zero during the simulation. The final configurations 
obtained are good examples of the typical morphologies 
favored by the competition of the different energy terms. 

Results will be presented for a unique set of values of 
the coefficients a and 7 in the gradient and substrate 
interaction terms [Eqs. Q and (Jj)]. It can be shown 
that a change of these parameters can be absorbed in 
a rescaling of the spatial coordinate and a global redef- 
inition of the free energy |7j|. We have chosen the values 
a = 3, 7 = 0.01, for which the expected spatial scale of 
the structures to be seen (based on the estimation in the 
previous section) is of the order of ten mesh parameters. 
Thus, the geometrical features of the patterns will be 
reasonably larger than the discretization of the numer- 
ical mesh, and this is appropriate to eliminate spurious 
effects associated to this discreteness. 

The configurations obtained are presented in Figs. [21 
to El corresponding to cases A-to-D, respectively. Each 
panel in each figure represents a different value of ei, 
i.e., a different degree of mismatch with the substrate. 
Each main plot shows the spatial distribution of e\ or 
ed (depending on which is more representative in each 
case). The scale for this plot is indicated in the bar of 
the inset. The inset also shows the combined distribution 
of ed vs. e\ for all elements in the system. We now discuss 
separately each case. 

Except for the existence of the substrate, case A is 
qualitatively analogous to a well studied case, which is 
applicable for instance to phase separation in alloys H,||. 
In the case in which the collapsed state has the same 
shear modulus than the original state (case Al), there 
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is strong evidence that the possible ground states corre- 
spond to a striped phase, a bubble phase, or a uniform 
phase depending on the value of e±. The patterns ob- 
tained in the numerical simulations presented here (Fig. 
I2J are in reasonable agreement with these results, al- 
though it is clear that the perfect striped phase is not 
easily obtained. The reason for this discrepancy is that 
there is a large gain in entropy when the stripes disorder, 
and it is very difficult numerically to get rid of this effect. 

In the case in which the collapsed phase has a lower 
shear modulus than the original phase (case A2, Fig. 
0) , bubbles of the softer phase tend to be unstable with 
respect to elongation^ [|J. This is obvious in Fig. [3] 
compared to Fig. Note that the effect is not sym- 
metric: if the mismatch is such that the minority phase 
corresponds to the more rigid phase, then the bubbles 
remain stable. In both cases Al and A2, for even larger 
mismatch (not shown) all the system becomes uniformly 
collapsed, with ei(r) = ei, e2(r) = es(r) = 0. 

For case B (Fig. 0J) we observe a stronger tendency 
to form stripes and polygonal patterns. There is also 
an important difference with case A2 for large mismatch 
(lower-right panel): Now each element prefers to be as 
close as possible to the collapsed minimum, that locates 
at a non-zero value of ed- However, not all elements can 
have the same values of ei and e^, since the spatial av- 
erages e-i and £3 should be zero. The existence of the 
substrate also discourages uniform phases with constant 
e2 or 63. The configuration of the system is such that in- 
dividual elements tend to be distributed close to the ring 
of minima in the ei-e-i plane. In real space, singularities 
appear at which e2 ~ — 0, and around them the con- 
figuration point in the e2-e3 plane rotates 2ir, clockwise, 
or counterclockwise, thus defining 'vortex-like' or 'anti- 
vortex-like' defects. Density of these defects is mainly 
controlled by the strength of the interaction to the sub- 
strate (density tends to zero for vanishing interaction). 

Cases C and D (Figs. EJandEJ) present an interesting 
difference in morphology with respect to previous cases: 
Stripes of collapsed regions present now a wavy struc- 
ture. In some cases, we observe even the existence or 
ring shaped collapsed regions. The origin of wiggling 
collapsed stripes is of course dictated by the tendency of 
the model to minimize the energy. An elemental demon- 
stration of the instability of straight collapsed regions 
with respect to wiggled patterns will be presented now. 
Consider the sketch of Fig. |7Ja) of a straight collapsed 
region along the x direction. The arrows in the plot are a 
schematic representation of the displacement field u(r). 
Suppose that we have determined the best form of the 
function u(r) = uo(y)y in order to minimize the energy. 
We will see that in some circumstances the energy can be 
reduced further by changing of u(r) = uo(y)y to a new 
displacement field of the form u(r) = uo{y — As'm(kx))y. 
Note that this represents a wiggling pattern with oscil- 
lation amplitude A and wave vector k (as sketched in 
Fig. db)). By calculating the form of our fundamental 
variables ei, e2 and e3, and plugging these values into 



the different term of our total energy, we see that the 
substrate energy does not change at all. The gradient 
term is increased in a quantity (per unit length in the 
x direction) of the order of A 2 k 2 d, where d is the thick- 
ness of the collapsed region. The local term changes only 
because of the existence of a non-zero value of e%. This 
change is also proportional to A 2 k 2 d, and the overall fac- 
tor can be negative (i.e., we can have an energy gain) if 
v (= e^/(ef — e")) is sufficiently large, as depicted in 
Fig. EIc). Then we see that if this quadratic decrease 
of energy overcomes the quadratic increase due to the 
gradient terms, a wavy stripe is expected. This lowest 
order sketchy argument does not provide values for the 
amplitude A and wave vector k of the distortion. They 
could in principle be obtained by considering higher or- 
der terms in the expansions. We failed in doing such a 
calculation. We only notice that from the results of the 
numerical simulations (Figs. |^|and|^)l it seems plausible 
that both A and the wave length k^ 1 are of the order of 
the stripe thickness d. 



IV. PRACTICAL REALIZATIONS 

We have analyzed an idealized situation of a material 
that has two well defined elastic configurations of mini- 
mum energy. We have shown that non-trivial morpholo- 
gies may appear when a quasi-two-dimensional piece of 
such a material is uniformly attached to a rigid substrate. 
In searching for practical realizations, it will be difficult 
to find that all the above assumptions are satisfied, how- 
ever, our model can be a good idealized case to study 
some practical problems. Two examples will be presented 
now. 

The first case corresponds to the patterns observed in 
delamination of thin films|inl|. In this case, typically, a 
film is grown onto a substrate, and due to chemical dif- 
ferences between substrate and film (which may be en- 
hanced by temperature changes) elastic stresses develop 
between them. If the stresses are compressive into the 
film, this may induce the buckling up of part of the film, 
giving rise to patterns that have been termed 'telephone- 
cord' like, due to their wavy appearance. In some cases, 
wavy polygonal patterns of buckled regions are observed, 
remarkably similar to those in our Figs. |5l and |HJ In de- 
lamination patterns the wave length and amplitude of 
the undulations observed are both of the order of the 
thickness of the buckled region. 

This delamination process has some features that are a 
priori very different from the assumptions in our model. 
Two things look in fact very different: in delamination, 
the three dimensional nature of the process plays an im- 
portant role. In contrast, our model considers a two di- 
mensional geometry. Secondly, in delamination the in- 
teraction to the substrate is lost in the buckled regions, 
whereas in our model this interaction is always present. 
However, the similar appearance of some of our results 
with the delamination phenomenology is very suggestive. 
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We think that a justification of the similarity can be 
based on the following considerations: a stripe of buck- 
led material is qualitatively depicted in Fig. E^a). K we 
want to make an effective two-dimensional representation 
of this buckled film, we should associate to the buckled 
region a piece of the two dimensional model with higher 
density. Note also that the orientation of the buckled 
stripe can be any in the x-y plane, but that the contrac- 
tion that induces the density increase occurs only per- 
pendicularly to the buckled ridge. This means that we 
can associate buckling with our collapse transition with 
v — \. In delamination, the possibility of buckling de- 
pends on the possibility of detaching from the substrate. 
Although in our model the interaction with the substrate 
is always present, the free energy with two minima takes 
phenomenologically into account the two configurations 
of the film: attached and unbuckled (represented by the 
isotropic minimum of our model) and detached and buck- 
led (represented by the finite collapsed state). 

A second problem to which our model can be applied 
is the epitaxial growth of competing phases with differ- 
ent crystallographic parameters^]]. I n fact, if the grow- 
ing is coherent (to some extent), elastic stresses will be 
accumulated during growth, and which phase is chosen 
locally to grow the material will be dictated by a ten- 
dency to minimize the total elastic energy of the system. 
Thus, our model with two minima must be thought as 
describing the two different phases that can be chosen 
in the growth process. But in this case one of the pre- 
viously used assumptions needs to be changed: epitaxy 
is strongly dependent on the crystalline structures of the 
phases, and isotropy of the model in the x-y plane is not 
justified. However, the model can be easily modified to 
account for anisotropics in the x-y plane. We will present 
here only a single example of the possible outcome ob- 
tained in an anisotropic model. Let us consider that the 
substrate has a square symmetry (assumed to be aligned 
with the x-y coordinate axis), that one of the growing 
phases has also a square symmetry with a larger unit 
cell volume, and that the second growing phase has a 
rectangular symmetry with a lower unit cell volume than 
that of the substrate. This situation corresponds to a 
local free energy /o having an elastic minimum as usual, 
and a collapsed minimum localized now at e2 = ±e^, 
63 = 0. In this way the ring of minima in the e2-e^ 
plane of the previous cases is replaced by a couple of 
minima. The exact form of /□ we use is given in the Ap- 
pendix. The morphologies obtained are shown in Fig. [5] 
As expected, the existence of a rectangular phase reflects 
clearly on the morphology of the patterns. Note how 
upon change of the fraction of rectangular phase (which 
is dictated, as usual, by the degree of mismatch with 
the substrate), the morphology changes from stripes of 
rectangular phase parallel to the axis (for low density of 
rectangular phase) , to stripes of rectangular phase along 
the diagonals (for the case of large fraction of rectangu- 
lar phase). It is worth to be noticed that the present 
case can also be interpreted in the context of a model 



recently introduced to describe a square-to-rectangular 
martensitic transformation^, 0. In fact, we can con- 
sider the isotropic minimum as representing the austen- 
ite phase, and the collapsed minima as the martensite 
phase, in its two different variants: — ie^. In the 
case of Refs. 0, (in which a substrate has not been 
included) , volume change during martensitic transforma- 
tion is assumed to be negligible, and a phase of diagonal 
stripes similar to that in the lower plots of Fig. [9] has 
been observed 0, [B| (note that in our case the width of 
the diagonal stripes is fixed by the strength of the inter- 
action with the substrate). From the results of the first 
plots of Fig. |5] we can conclude that a phase of stripes 
along the x and y directions is favored instead when the 
transformation has an appreciable volume change and a 
finite fraction of the austenite phase remains in the sys- 
tem. 



V. SUMMARY AND CONCLUSIONS 

Expansion ridges are analog to surface cracks, with the 
difference that they appear when a film expands with 
respect to the substrate (contrary to cracks, that ap- 
pear when the film contracts with respect to the under- 
lying material) . A non- linear elastic behavior of the film, 
namely the possibility of a collapse upon applied stress is 
necessary for expansion ridges to appear. We have sim- 
ulated expansion ridges through an elastic model that 
uses the components of the strain tensor as fundamen- 
tal variables. We have shown how complex morphologies 
can appear due to competitive elastic interactions in the 
system. By changing the parameters of the model, col- 
lapsed regions in the form of bubbles, straight stripes, un- 
dulating stripes and rings have been obtained. We have 
argued that the model can be applied to understand the 
characteristics of delamination patterns in thin films, and 
bi-phase epitaxial growth. We have concentrated in the 
description of the morphologies, and then our presenta- 
tion has been qualitative to a large extent. A few impor- 
tant things can be highlighted: The relevant spatial scale 
of the patterns has been seen to be fixed by a competi- 
tion between interaction with the substrate and gradient 
effects. In the case in which the collapse is isotropic, 
we have qualitatively reproduced well known results in 
the field of alloy decomposition. Wiggle stripes of col- 
lapsed material appear when the collapse of an elemental 
circular piece of material is such that contraction in one 
direction is accompanied by expansion in the perpendicu- 
lar direction. We have compared these wiggling patterns 
with those observed in delamination of thin films. We 
have also provided an example showing that anisotropy 
can be easily incorporated in the model, and showed that 
the anisotropic case may have relevance for the study 
of bi-phase epitaxial growth. More detailed studies an 
comparisons on each of these individual realizations is in 
progress. 



5 



VI. APPENDIX 



We give here the explicit expressions for the local free energy f Q used in each of the simulations presented in the 
paper. For Figs. El and El we use 



fa 



i + (6i - I) 2 



[1 + tanh(3ei 



Cx^ + (e 1 + l) 2 



[1 — tanh(— 3ei) 



with d = 1 in Fig. El and C x = in Fig. El In Figs. H El El we use 



fo 



f + f + (ei - C 3 ) 2 



[l + tanh(C 2 ei)] + 



-y + f + (ei + C 3 ) 2 + i 



[1 - tanh(-C 2 ei)] 



(9) 



(10) 



with C 2 = 2 and C 3 = 0.4 in Fig. |C 2 = 3.5 and C 3 = 0.2 in Fig. El and C* 2 = 4 and C 3 = 0.1 in Fig. El Finally, in 
Fig. El we use 



P 2 4 
2 , e 2 , / ^ \2 

Y + -j- + (ei - C 3 ) 



/o = 

with C 2 = 2 and C 3 = 0.4 



[l + tanh(C 2 ei)] 



G 2 i e 2 , /„ , /^i \2 , 1 



-f + f + (ei + C 3 



[1 -tanh(-C 2 ei)] +e 3 



(11) 
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FIG. 1: Different possibilities for the location of the collapsed 
minimum in the ei-e^ plane. The elastic minimum is located 
at 0. The sketches show qualitative an originally circular piece 
of material at O and its collapsed form for different locations 
of the collapsed minimum. In cases A and B the an original 
circular piece of material contracts in all directions, whereas 
for C and D there is a direction in which it actually expands 
(chosen here to be the vertical direction). 
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FIG. 2: (Color online) Spatial configurations obtained for case 
Al, for values of ei equal to —0.4, —0.2, —0.1, and 0.2 (from 
left to right and from top to bottom). The shadows indicate 
values of ei, with the scale indicated in the inset. In the inset 
plots, the distribution of ea vs. ei is shown for all elements 
of the system. 




FIG. 3: (Color online) Same as Fig. EJfor case A2, and values 
of ei of -0.4, -0.2, -0.1, and 0. 
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FIG. 4: (Color online) Same as Fig. [5]for case B, and values 
of ei of —0.2, 0, 0.4, and 0.5 (note that in this and next two 
Figures, the main plot shows values of e<j, instead of ei). 




FIG. 5: (Color online) Same as Fig. 0]for case C, and values 
of ei of -0.05, 0.1, 0.3, and 0.4. 



FIG. 6: (Color online) Same as Fig. 2]for case D, and values 
of ei of -0.005, 0.05, 0.2, and 0.4. 
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FIG. 7: (a) The structure of a straight collapsed region. Ar- 
rows indicate qualitatively the displacement u(r). This con- 
figuration can be unstable with respect to the deformation 
shown in (b), for certain forms of the local free energy /o. In 
fact, part (c) illustrates how a piece of material collapsed in 
a single direction (represented by the dotted circle) can lower 
its energy by generating a finite shear distortion es if this 
brings it closer to the ring of minima of the potential. 
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FIG. 8: A stripe of delaminated thin film, and the qualitative 
equivalent as a region of collapsed material in our model. 




FIG. 9: (Color online) Morphologies for an anisotropic case in 
which the collapsed minima are located at e2 = ±e§ ', e-z = 0. 
The values of e\ for the plots are —0.1, 0, 0.1, 0.3, 0.4, and 0.5 
respectively. Note that as long as there is some fraction of the 
system in the isotropic elastic minimum (first three plots) the 
morphology of the collapsed regions is that of inter-crossed 
stripes along x and y, while when the fraction of the system 
in the elastic configuration is negligible (last two plots) the 
morphology changes to stripes along the diagonals. 



